Theoretical analysis of Cooper-pair phase fluctuations in underdoped cuprates: 

a spin-fluctuation exchange study 



On 

o 

Q 



o 
o 

di- 
ss 

S 

i 

-a 

£ 

o 
o 



(N 
> 

\o 

O 

(N 

r-H 

^— > 



i 

c 

o 
o 



X 



D. Manske 1 , T. Dahm 2 tl, and K.H. Bennemann 1 

1 Institut fur Theoretische Physik, Freie Universitat Berlin, Arnvmallee 14, D-14.195 Berlin, Germany 
2 Max- Planck- Institut fur Physik komplexer Systeme, Nothnitzer Str. 38, D-01187 Dresden, Germany 

(February 1, 2008) 

We study Cooper-pair phase fluctuations in cuprate superconductors for a spin fluctuation pairing 
interaction. Using an electronic theory we calculate in particular for the underdoped cuprate su- 
perconductors the superfluid density n s (T), the superconducting transition temperature T c (x) oc n s 
below which phase coherent Cooper-pairs occur, and T* (x) where the phase incoherent Cooper-pairs 
disappear. Also we present results for the penetration depth X(x,T) and for the weak pseudogap 
temperature T*(x) at which a gap structure occurs in the spectral density. A MeiBner effect is 
obtained only for T < T c . We find that phase fluctuations become increasingly important in the 
underdoped regime and lead to a reduction of T c in good agreement with the experimental situation. 
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In accordance with the experimental observation for 
underdoped cuprate superconductors that the supercon- 
ducting transition temperature is related to the super- 
fluid density n s and behaves as T c oc n s (0) , in contrast to 
the overdoped case |ffl] , we use an electronic theory which 
takes into account phase fluctuations of the Cooper-pairs. 
From this one expects a phase diagram with T* at which 
phase incoherent Cooper-pairs occur and a lower transi- 
tion temperature T c oc n s at which these Cooper-pairs 
become phase coherent. Such a physical picture for the 
underdoped cuprates has been treated phenomenologi- 
cally before by Chakraverty et al. 0, more recently by 
Emery and Kivelson || , and also by Schmalian et al. [Q . 

Close to T c one gets as in Ginzburg-Landau (GL) the- 
ory the free energy (without external field) 
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where AF BCS = ±N(E F )A 2 and AF phase = ^n s with 
n s = n s (V$(r)V < l'(0)). Here, we assume only a spatial 
dependence of the phase <&(r) of the superconducting or- 
der parameter. This implies a superconducting transition 
at k B T c = AFbcs, when AF B cs < AF phase and Coper- 
pairs are formed, and ksT c — AF p hase cx n s (0), when 
AFphase < AFbcs and Cooper-pairs become phase co- 
herent. This neglects 2d effects. In this letter we use 
an electronic theory in order to calculate the doping de- 
pendence of N(Ep), Ao(a;), and n s (x) and thus the two 
contributions to AF and the two doping regimes for non- 
BCS and BCS behavior of the superconducting transi- 
tion temperature. Based on a spin fluctuation exchange 
mechanism we determine the doping dependence of the 
phase diagram with T*(x), T c (x), and also T* at which 
far above T c a gap appears in the spectral density. In 
accordance with Eq. (|j) we find from our microscopic 
calculation that phase fluctuations of the Cooper-pairs 
are significant for the underdoped cuprates. For the over- 
doped superconductors we obtain a mean-field behavior 
and T^ xp ~ T*. Our results suggest (see also Kosterlitz- 
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FIG. 1. Phase diagram for high-T c superconductors result- 
ing from a spin fluctuation induced Cooper-pairing includ- 
ing their phase fluctuations. T c * denotes the temperature be- 
low which Cooper-pairs are formed. The dashed curve gives 
T c oc n s (T = 0,:r). Below T* (see inset) we get a gap struc- 
ture in the spectral density. 

Thouless (KT) theory) that a MeiBner effect occurs only 
for T < T c (x), since A -2 oc n s . 

We treat the cuprates by using as a model the 2D Hub- 
bard Hamiltonian for a Cu02-planc 



H 



U Vn iT n u . (2) 



Here, creates an electron with spin a on site i, U 
denotes the on-site Coulomb interaction, and ty is the 
hopping integral. From H and the functional differenti- 
ation of the free energy F with respect to the Green's 
function Q, 6F{H}/SG = S, one obtains with the help of 
the Dyson equation the self-energies. Note, this gives in 
general coupled equations for the amplitude and phase 
<f>(r) of Q H and in accordance with the GL (Wannier) 
type treatment an energy gain due to phase coherence 
of Cooper-pairs and due to Cooper-pair condensation 
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It follows from Eq. (g) that the quasiparticle self- 
energy components X v (y — 0, 3, 1) with respect to 
the Pauli matrices t v in the Nambu representation |t]] 
are given by 

/>oo 

X u (k,oj) =N~ 1 V / rffi [P s (k-k',Q)±P c (k-k',fi)] 
t/ Jo 



x / &t//(u>,Q,a/)A,(k',c/) 
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Here, we used the T-matrix | 
change (FLEX) approximation 
Bcrk-Schrieffer-like 13) pairing 
structed with the dressed one-electron Green's func- 
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fluctuation interaction is given 
U 2 Im (3xs - Xso) with Xs = 
and the charge fluctuation interac- 
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= (27T)- 1 

uxsor 1 

Pc = (2vr) 1 U 2 Im (3 Xc - Xco) with X 
-l 



UXco) \ where Im x s o,co(q, ^) is given in Ref. 
|l2| . The subtracted terms in P s and P c remove a double 
counting that occurs in second order. In Eq. (|3|) the plus 
sign holds for Xq (quasiparticle renormalization) and X% 
(energy shift), and the minus sign for X\ (gap param- 
eter). The kernel / and the spectral functions A v are 
given in Ref. 

For the numerical evaluation we use a bare 
tight-binding dispersion relation e(k) = 2t [2— 
cos(7rfc x ) — cos(7rfc a ) — /i] and U/t — 4. Then, the dop- 
ing dependence n = ^ k = 1 — x is determined with 
the help of the k-dependent occupation number — 
2 f_ du!f(ui)N(k,uj) that is calculated self-consistcntly, 
where N(k,uj) — Ao(k, uj) + ^(k, u>). n = 1 corresponds 
to half filling. Our numerical calculations are performed 
on a square lattice with 256x256 points in k space of the 
Brillouin zone and with 200 points on the real u) axis 
up to 16t with an almost logarithmic mesh. The full 
momentum and frequency dependence of the quantities 
is kept |14|. Thus, our calculation includes pair break- 
ing effects on the Cooper-pairs resulting from lifetime 
effects of the elementary excitations. T* is determined 
from the linearized gap equation and the superconduct- 
ing state is found to have d-wave symmetry fL2| . The 
transition temperature T c at which phase coherence oc- 
curs is determined by T c oc n s (x) where the superfluid 
density n s (x,T)/m is calculated self-consistently from 
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Ss is the corresponding expression in the superconduct- 
ing state. Here, we utilize the f-sum rule for the real part 
of the conductivity CTi(w), i.e., / °° cri(w) duj = ire 2 n/2m 
where n is the 3D electron density and m denotes the ef- 
fective band mass for the tight-binding band considered. 
<r(w) is calculated in the normal and superconducting 
state using the standard Kubo formula |Tq |. Vertex cor- 
rections have been neglected. Physically speaking, we 



are looking for the loss of spectral weight of the Drude 
peak at to — that corresponds to excited quasiparticles 
above the superconducting condensate for temperatures 
T < T* . Finally, the penetration depth X(x,T) is calcu- 
lated within the London theory through n s cx A~ 2 ]l6| ]. 

In Fig. ^] we show the phase diagram for the cuprate 
superconductors obtained from our electronic theory and 
compare with the generalized experimental phase di- 
agram T c (x) that describes many superconductors as 
pointed out by Tallon Q. We get T c w T* for 
x > x opt , whereas for the underdoped superconductors 
T c oc n s (T — 0) agrees much better with the experimen- 
tal results than T* . Note, for the overdoped cuprates 
one gets the superconducting transition at T*, where 
phase coherent Cooper-pairing occurs since AF p h aS e be- 
comes largest. T* results from 4>{x,T*) = 0, where 
4>{k,ui) — Z(k,uj)A(k,Lj) is the strong-coupling super- 
conducting order parameter and is a mean-field result 
in the sense that Cooper-pair phase fluctuations have 
been neglected. However, the AF fluctuations are treated 
well beyond the mean-field approximation. It is remark- 
able that we find approximately that T c varies linearly 
with x. This cannot be the case for n s (x — * 0). We 
also expect that an improved inclusion of the pseudogap 
behavior and reduction of the density of states (DOS) 
due to antiferromagnetic (AF) correlations will yield a 
faster decrease of T c (x) for x — > 0. Also important is 
the fact that we find an optimal T c for x ~ 0.15 since 
for x < x op t we find T c oc n s (T = 0) and for x > x op t 
one has T c oc A(T = 0). Of interest are also the results 
for the pseudogap temperature T* at which a gap struc- 
ture in the DOS in the normal state occurs. This has 
been discussed in Refs. [P^fl- Such a behavior is seen 
experimentally in various SIN tunneling data, reflectiv- 
ity measurements as well as in the two-magnon response 
in Raman scattering Jl^,[l{| . Physically speaking, below 
T* the electrons are strongly coupled to the magnetic de- 
grees of freedom which are generated due to paramagnon 
(spin fluctuation) excitations in the system. Finally, we 
would like to mention that at x* = 0.19, T* and T* co- 
incide. Furthermore, our results indicate that for x > x* 
the non-BCS behavior T c oc n s (T = 0) seems not valid 
anymore. In contrast to the underdoped regime, there 
the energy gain due to Cooper-pair formation is smaller 
than the energy to break up phase coherence. Thus, 
Cooper-pair phase fluctuations are unimportant in this 
regime. Note, T c may be viewed in analogy to the Curie 
temperature in Ferromagnets like Ni or Fe, where spin- 
disorder occurs. 

In order to investigate different models used in the con- 
text of Cooper-pair phase fluctuations J^-|J|, we discuss 
now the superfluid density for finite temperatures. In 
Fig. |^(a) results are given for n s {T,x)/m. To obtain 
n s (T = 0) we have used a polynomial fit up to third 
order ^0|. For the doping values investigated here, we 
find a linear behavior of n s (0,x) except for the strongly 
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FIG. 2. (a) Results for the superfluid density divided by m 
in units 1Q~ 2 m'^ 1 A~ 2 as a function of temperature T for var- 
ious doping concentrations: x=0.221 (open circles), x=0.201 
(circles), x=0.188 (triangles), x=0.155 (squares), x=0.122 (di- 
amonds), x=0.091 (circles). n a vanishes at T = T c *. From 
n 3 (x,T), T c is given as indicated by the intersection with 
aT, where 1/a = 2.202 (3D XY), 1/a = n/2 (KT theory), 
1/a = 0.9 (2D XY). The obtained results for T c are shown in 
(b). 

overdoped case, i.e. x > 0.22 where n s (T = 0) starts 
to decrease. This behavior indicates a crossover to the 
(BCS) weak-coupling limit where n s = n for T = is ex- 
pected. This implies that A(T = 0, x) is asymmetric with 
respect to optimal doping as suggested from Fig. [|. Note 
that the energy AF p hase oc n s (0) is still larger than the 
energy gain due to Cooper-pair condensation. Thus we 
conclude that Cooper-pair phase fluctuations are unim- 
portant in the overdoped regime. Note, n s (T,x) — > for 
T — > T*, since Cooper-pairs disappear at T*. However, 
the phase coherence temperature T c has to be determined 
by spatially averaging over the Cooper-pair phase fluctu- 
ations. In the presence of spatial phase fluctuations the 
average superfluid density n s will vanish at T c so that as 
in KT theory no Meifiner effect occurs above T c . Note, 
within KT theory T c is given by 

1 % 2 n s{T c ) . . 

- n — = k B T c , (4) 

a 4m 



where 1/a = w/2. In the case of the similar 2D XY 
model one has 1/a = 0.9. Our construction used in Fig. 
^ involves only the determination of n s (T = T c ), but does 
not imply a jump in n s at T c which is true only for the 
2D XY-mode\ and KT theory. For a comparison with 
Ref. [Q, we show also results for the 3D version of the 
XY-modcl. All these theories take into account phase 
and amplitude of the Cooper-pairs. We see that in the 
2D (3D) version of the XY-mode\ T c is lower (higher) 
compared to the KT theory. Physically speaking, in the 
AY-model n s (or the phase stiffness) is the only rele- 
vant energy scale. Thus one always has T c oc n s . In 
the KT theory one has two energy scales, namely the 
phase stiffness and also the vortex core energy. Never- 
theless, all these theories yield within our microscopic 
treatment an optimal doping concentration. Note, we 
obtain n s (x — > 0) due to pair breaking resulting from the 
decreasing lifetimes of the Cooper-pairs |2E|. Improved 
calculations taking into account AF correlations leading 
to static antifcrromagnetism cause a further decrease of 
n s {x) for i-f0 |§. 

In Fig. |^ we show the mean field results for the pen- 
etration depth. In (a) the overdoped case is displayed: 
our results for X~ 2 (T) agree qualitatively with the ex- 
perimental data [^2],^]]. It is remarkable that we find 
a linear behavior of A _3 (T) (see inset). In Ref. 0] at 
T = 0.9T* a value for A 3 (0)/A 3 (T) of 0.2 has been found. 
Instead, we find 0.16. Note that we did not include crit- 
ical fluctuations. A closer inspection of our data leads 
to the conclusion that the rapid opening of the super- 
conducting gap -which is calculated self-consistently- is 
the main reason for the reported behavior of the penetra- 
tion depth (T^j2^]. Thus, we conclude that the existence 
of critical fluctuations is not necessary in order to un- 
derstand the observed behavior of the penetration depth 
close to the critical temperature in the cuprates. This 
might shed some light on the significance of vertex cor- 
rections. However, very close to T c we cannot present 
numerical results yet. For this range it remains unclear 
whether or not a change of the powerlaws from 1/3 to 1/2 
occur. Fig. |](b) indicates that for calculating A(T, x) one 
must use the superfluid density n s referring to phase co- 
herence as is done in KT theory. Note, static KT theory 
predicts an universal jump in n s at T = T c as indicated, 
which is not observed experimentally j2^-||J]; the n s re- 
sults obtained from our electronic theory refer only to the 
Cooper-pair density n s (x, T) but n s — n s at T < T c . As 
mentioned earlier, the spatially phase averaged n s (x,T) 
yields n s (T c ,x) — for incoherent Cooper-pairs. The 
comparison with experiment indicates that n s < n s al- 
ready for T < T c . No Meifiner effect is expected for 
T c < T < T*. While our calculation of n s (T, x) is ex- 
pected to give the correct tendency for the behavior of the 
underdoped superconductors, stronger AF effects must 
be included for x — > 0. Also, quantum phase fluctuations 
are neglected. All this would cause a more rapid decrease 



3 



(N 



1.0 




0.9 




0.8 




0.7 




0.6 


0.20 


0.5 


0.16 


0.4 


0.12 


0.3 


0.08 
0.04 


0.2 


O.Ofl 


0.1 




0.0 






0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

t/t: 



r< 

o 

IN 



0.0 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

T / Tc 

FIG. 3. Results for the penetration depth X(T, x) 
(n B oc A -2 ). The inset in (a) shows a remarkable linear be- 
havior of A~^(T) for T ~ T* even without having included 
critical fluctuations. In (b) we indicate the behavior expected 
for n s and derived from KT theory, where n s (T) — > at T c . 



of T c oc n s (x) in Fig.|l|. 

In summary, it is interesting that our electronic theory 
using the model Hamiltonian H, Eq. (^), gives remark- 
able agreement with some basic experimental observa- 
tions. In particular, we have shown within a microscopic 
theory based on a spin fluctuation pairing mechanism 
that Cooper-pair phase fluctuations become important 
for underdoped cuprates. This gives a microscopic justi- 
fication for the phenomenological theories in Refs. [§,1). 
Of course, a fully self-consistent determination includ- 
ing the phase fluctuations of the Green's functions and 
the electronic exitations should be performed | p6[ . Thus, 
one might get also an insight about the origin of the 
pseudogaps above T c and revealing asymmetric doping 
dependence of various properties with respect to x op t- 
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